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Abstract 



The paper presents an analytical description of the growth of a two-component bubble 
in a binary liquid-gas solution. We obtain asymptotic self-similar time dependence of the 
bubble radius and analytical expressions for the non-steady profiles of dissolved gases 
around the bubble. We show that the necessary condition for the self-similar regime of 
bubble growth is the constant, steady-state composition of the bubble. The equation for 
^ I the steady-state composition is obtained. We reveal the dependence of the steady-state 

composition on the solubility laws of the bubble components. Besides, the universal, 
independent from the solubility laws, expressions for the steady-state composition are 
obtained for the case of strong supersaturations, which are typical for the homogeneous 
nucleation of a bubble. 



> 

Introduction 

o 



The subject of this paper is theoretical treatment of growth via diffusion of a two- component 
gas bubble in a supersaturated gas-liquid solution. Understanding the regularities of gas bubble 
growth in gas-liquid solutions is crucial for the study of various natural phenomena and for 
the management of different technological processes. Gas bubble growth is important e. g. 
for glass refining processes [1], polymeric foams production [2]. The growth of gas bubbles 
\ dissolved in magma is a process which governs volcanic eruptions [3-6]. 

The problem of diffusional bubble growth with single gas in it was first described theoreti- 
cally in the classical paper by Epstein and Plesset [7] and remains topical till nowadays when 
treated under different specific conditions [8-14]. 

While the discussion in Refs. [1,2,5-14] is limited to bubbles containing only one gas, often 
there are two or more different gases dissolved in the solution and, therefore, presenting in the 
bubble. The growth of a multicomponent bubble in glass melts was described theoretically 
by Ramos [15] and Cable and Frade [16]. The main question in these studies was the time 
dependence of multicomponent bubble radius. In Ref. [15] the growth was considered from a 
small initial size (when Laplace pressure in the bubble is of the order of the external pressure 
of the solution or more, and therefore, the surface tension strongly influences bubble growth); 
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and there the problem was treated numerically only. However, Ramos [15] reported a small 
effect of surface tension on bubble growth. 

Cable and Frade [16] neglected the effect of the surface tension initially. They solved the 
problem numerically for the arbitrary (but large enough to neglect Laplace pressure) initial 
radius and also presented an analytical description of the growth from the zero initial radius. 
Ref. [16] states that the composition of a multicomponent bubble tends to a steady-state value; 
however, explicit equations for calculating this composition have not been obtained there. (In 
Ref. [16] this value of the gas composition in the bubble was erroneously called "equilibrium". 
If the composition was equilibrium, the bubble would be unable to grow). Besides that, in 
Refs. [15, 16] the discussion was limited to low supersaturations typical for problems related 
with glass production. Finally, the influence of gas solubility laws on the steady composition 
of the bubble was totally excluded from discussion. 

In the current paper we will consider a bubble which consists of two different gases in a liquid 
solution where both of these gases are dissolved. Still, some of the results can be generalized 
for the case when the number of components exceeds two. We will obtain an asymptotic 
self-similar time dependence of a bubble radius and analytical expressions for the non-steady 
profiles of dissolved gases. We will also obtain the equation for the steady-state composition 
of the bubble. The steadiness of the composition is the necessary condition for the self-similar 
regime of bubble growth. The dependence of the steady-state composition on the solubility 
laws of the components will be demonstrated. For the case of strong supersaturations (which 
is typical for the homogeneous nucleation of a bubble) we will obtain universal expressions 
for the steady-state composition that are independent from the solubility laws. Since the self- 
similar approach to the bubble growth description will be used, we will consider relatively large 
bubbles. The self-similar regime of bubble growth takes place only when the influence of the 
surface tension could be neglected. 

The diffusional growth of a two-component bubble is similar to the growth of a two- 
component liquid droplet in a supersaturated vapor-gas mixture [17-19]. Nevertheless, there is 
a significant difference between these two processes. While homogeneous nucleation of droplets 
in supersaturated vapor occurs at relatively low supersaturations (of the order of a unit) and 
vapor diffusion toward a growing drop can be considered steady [20], nucleation of gas bubbles 
in supersaturated solutions usually takes place at high values of supersaturation; therefore, as 
the radius of the bubble increases, the steady regime of diffusion gradually gives way to the 
non-steady one [13,14]. As we will show further, the steady-state composition of the bub- 
ble depends on the degree of steadiness of diffusion fluxes of dissolved gases molecules to the 
growing bubble. 

2 Self-similar diffusion problems 

The state of the solution is stipulated by temperature T, pressure II and the initial densities 
711,0 and 77-2,0 of the dissolved gases (hereinafter we will refer to number density of molecules 
using the term "density"). A gas bubble nucleates in the solution and starts growing regularly. 
We will consider the bubble after some time from the moment of its nucleation, i. e. when its 
radius R has reached such a value that it satisfies strong inequality 

R > 2<r/n. (1) 

Quantity a here is the surface tension of the pure liquid solvent (we consider only diluted 
solutions without dissolved surfactant). Eq. (jTJ) means that the influence of Laplace forces on 
the bubble growth dynamics is negligible; and the pressure in it is equal to the pressure of the 
solution II. 
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The composition of the bubble is determined by concentrations of both gases in it. Under 
the notion of concentration we understand the mole fraction: 



Ni 



a = 



^ + N 2 



Hereinafter, when index i is used, we mean that the value corresponds to both cases i — 1 and 
i = 2; Ni and N 2 are the numbers of molecules of gases 1 and 2 in the bubble. Evidently Eq. 
O) is equivalent to 



n 



1 41 



ni,g + n 2: g 

where ni jS , n 2g are the densities of gases 1 and 2 in the bubble, 



(3) 



n i,9 



N,\fR* 



(4) 



When strong inequality ([I]) is fulfilled, the density profiles can be obtained by solving the 
following diffusion problems (see e. g. Refs. [12,21] ): 



drij(r,t) _ Ehd_ 



Of 



2 drii(r,t) 



Or 



dr 



n>i,o, 



(5) 



(6) 



ni(r,t)\ r=R = n i:00 (ci). (7) 

Here r is the distance from the bubble center; rii(r,t) is the 2-th gas density profile; Di is 
the diffusion coefficient of the dissolved i-th gas in the solvent (diffusion coefficients can be 
assumed constant for the diluted solution); R = dR/dt is the rate of change of the bubble 
radius in time. Quantity rij j00 (ci) is the equilibrium density of the i-th gas at the flat surface 
of the gas phase with the concentration Cj of the 2-th gas in it. 

Even in the form (JSD, ©, O the problem is non-trivial for analytical treatment. Neverthe- 
less, if we formally assume that Eqs. ([5]), (J6]), ([7]) are valid from t = and set the initial bubble 
radius value R\ t=0 to as in Ref. [16] , we can write self-similar solutions for diffusion problems, 
as it was previously done for the one-component case [12,21]. It has to be noted that after the 
pioneering works by Zener [22] and Frank [23] , who used self-similar method for the description 
of diffusional crystal growth in the supersaturated solution, this method was exploited repeat- 
edly by various authors for droplet growth in supersaturated vapor-gas medium [18,20,24] and 
for bubble growth in superheated and supersaturated solutions [8-12,21]. 

Following Refs. [24] and [12] , to find a self-similar solution we introduce dimensionless 
variable p as 

p = r/R. (8) 

and presume that rii(r,t) = rii(p). Then, Eq. for rii(p) can be obtained from Eq. §5§ and has 
the form 



d 2 rij(p) 
dp 2 



2 RR 



i 

p 2 



drii(p) 
dp 



0. 



(9) 



Initial and boundary conditions ([6]) and ([7]) also can be rewritten in terms of dimensionless 
variable p as follows 

rii(p)\ p=00 =n i>0 , (10) 



Mp)\ 



p =i 



IF 
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The necessary condition for Eq. (Q to have a self-similar solution is independence of fraction 
RR/Di from time, i. e. 

RR = const. (12) 

To assure whether the fulfillment of Eq. (I12p is possible, we need to consider the equation of 
material balance: 



dt dr 



(13) 



r=R 



Substituting Eq. (TjJ into Eq. (TTB"]) and transforming derivatives on r to derivatives on p [using 
Eq. USD], w e obtain 

1 drii(p) 



RR = Di 

n itg dp 



R 2 1 dn it9 ^ 



p=i 3 fii g dt 



In order for Eq. (fl4"|) to be in concordance with Eq. (fl2"l) . it needs to have the second addend 
in its r. h. s. equal to 0. It means that nj i9 is constant in time; and this situation corresponds 
to the steady-state composition of the bubble. We will denote the values of these constant 
densities (and constant concentrations Cj) with superscript "s" - steady-state: nf g and cf. 
Therefore, Eq. ( 1141) can be rewritten as 

RR = DA, (15) 
where the dimensionless parameter bi is introduced via 

1 drii(p) 



bi - 

n*, g d P 



(16) 

P =i 

Integrating Eq. Qj, using Eqs. f|T6|) and ffTTj) . we have 

The value of parameter bi in Eq. (1T71) is not known yet. To obtain this value we will use 
condition ([TO]) , which was not exploited earlier. Substituting Eq. ([TO]) in Eq. f[T7]) . we have 

6 . e 3»*/2 ^ e -M 2 /2-^/x_ (lg) 



where important dimensionless parameter is introduced using 



(19) 



Transcendental equation (TT5]) has the following asymptotics (see Ref. [12] ): 

k ~ a< (aj /2 < 1) (20) 

and 

6i - -a- (a, > 10). (21) 

71 

Eq. ([20]) corresponds to the situation when gas diffusion to the bubble is steady, which can 
take place only at low supersaturation (it is possible only when the bubble nucleates hetero- 
geneously). Eq. (|2"T]) corresponds to the strongly non-steady-state diffusion, which is typical 
for homogeneous nucleation of gas bubbles in supersaturated solution. For details considering 
steady-state conditions of gas bubble growth see Refs. [13,14] . 
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Concluding the current section we should note that describing the bubble growth by diffu- 
sion equation (jSJ) we made two simplifications. First, we neglected the enthalpy of gas dissolu- 
tion, assuming that the temperature of the bubble is equal to the temperature of the solution. 
The analysis of this assumption was made in Ref. [12] for the growth of one-component gas 
bubble, and it was shown that this effect is negligible even at relatively high supersaturation, 
when a,i < 20. 

Another simplification is related with the presupposed condition of mechanical equilibrium 
between the bubble and solution, in particular, we neglected the solvent viscosity. Its influence 
on bubble dynamics can be estimated using Rayleigh-Plesset equation (see e. g. Ref. [25] ). 
For the gas pressure in the bubble Pr we will have: 

2a R 

Pr = IL+ — +477-, (22) 

where rj is the dynamic viscosity of the solvent. The inertial terms in Rayleigh-Plesset equation 
are negligible for any reasonable bubble growth rate. In our examination we already neglected 
the second term in Eq. (1221) due to strong inequality ([I]). To neglect the third term we need 
the following strong inequality to be fulfilled: 

V « (23) 
Taking into account Eq. (|15[) . we can rewrite this inequality as 

n R 2 

"«4AV (24) 

Let us estimate the value in the r. h. s. of inequality (1241) . First of all, we need evaluation 
for the bubble radius R from strong inequality ([T]). Values of surface tension both for water 
and for volcanic systems [26] are a ~ 10 _1 N m -1 ; pressure II ~ 10 5 Pa, thus we have the 
minimal radius of the bubble i?~2x 10 -5 m. The value of 6j is determined by the value of 
cii and for Oj ~ 20 reaches fej ~ 800 (see Eq. (|2ip ). Taking 10 -10 m 2 s -1 as the estimation for 
Di, we can rewrite Eq. (pHj) as 

r] < 10 3 Pa s. (25) 

It should be noted that "common" liquids at normal conditions always satisfy this condition: 
for water we have t] ~ 10~ 3 Pa s and for glycerol t] ~ 1 Pa s [27]. Even for volcanic systems 
the values of viscosity that satisfy strong inequality fl25|) are quite typical when Si02 content is 
not too high (basalt, andesite and dactite melts) [3]. However, for rhyolite melts (~ 70% S1O2) 
viscosity can reach the values of 10 7 Pa s [26]; and therefore for rhyolite melts the effect of 
solvent viscosity has to be taken into account. 

3 Steady-state composition of the bubble 

From Eq. (fT5]) one can derive the parabolic relation between the radius and time: 

R = (2Dibit) 1/2 . (26) 

It has to be noted that numerical solutions of such diffusion problems showed that the assump- 
tion of parabolic relation between radius and time Eq. (1261) is reasonable as an asymptotic 
behavior for multi-component bubble for arbitrary number of gases in it [16]. 
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The equality of the r. h. s. of Eq. ([15]) for both cases i — 1 and i = 2 gives us 

£>!&! = £> 2 6 2 . (27) 

Eq. (]27p is the necessary condition for the self-similarity of the two-component bubble growth. 
And, as we have seen, the growth of the bubble can be self-similar only when its composition (ci 
and C2) is constant. Eq. f[2"T|) together with Eq. ffTSl) determine the steady-state composition of 
a bubble. Analytical solution of these eqations for different particular cases will be presented 
further in this section. 

It should be noted that the condition on steady-state composition of growing compound 
in its general form analogous to Eq. (1271) was obtained for the case of crystal growth in the 
supersaturated solution [23] and for the case of droplet growth in supersaturated vapor-gas 
medium [18]. 



3.1 Arbitrary supersaturations 

We will consider two characteristic cases: when both diffusion fluxes to the bubble can be 
considered steady, and when both diffusion fluxes are significantly non-steady. Let us begin 
with the case when gas diffusion can be considered steady, and we can use Eq. ( 1201) for both 
i = 1 and i = 2. Thus, from Eq. (I2T|) we have: 

D x a x = D 2 a 2 . (28) 

Substituting definition (fl9|) into Eq. (128|) . taking into account n ij9 = nf and Eq. ([3]) (which 
gives us C2 = 1 — c{), we obtain 

Pi 1 ~ cf = n 2fi - n 2 ,oo(l ~ cf) 
D 2 c\ m )0 - ni l0O (cf) 

Equation analogous to Eq. ( 129]) was obtained by Kulmala et al. [17] for the steady-state 
growth of a two-component droplet. Eq. ( 1291) gives us explicit values for the steady-state 
concentrations in the bubble c\ and 1 — c{ , which are unambiguously defined when the solubility 
laws for each gas are set. Solubility laws give us the explicit dependence of equilibrium densities 
n i,oo(ci) and 7^2,00(1 ~ c'l) on corresponding gas concentrations, e. g. Henry's law (which 
describes gases where molecules do not dissociate during dissolution) 

(30) 

and Sievert's law (which describes gases where molecules dissociate into two parts during 
dissolution) 

Here n ijOQ = n i oo (l) is the density of the i-th dissolved gas which is in the equilibrium with 
the bubble of the pure i-th gas. Substituting relations (1301) or (131!) into Eq. (129|) . it is possible 
to obtain the analytical solution of the latter. 

Evidently, the steady-state composition is different for different solubility laws. To demon- 
strate the difference between the values of steady-state concentration c\ for Henry's law and 
Sievert's law, in Fig. [T]we plot the value {c s 1Henry — cl Sievert ) /c s 1HenTy as a function of rii i0 /ni iOO 
for different values of n 2t0 /n 2tOO (for the sake of simplicity we put D\ = D 2 and ni >00 = n 2 ,oo)- 
One can see that the smaller supersaturations of both components are, the higher the difference 
between steady-state concentrations for different solubility laws is. In the end of this section 
we will demonstrate that when the solution is strongly supersaturated with both components, 



6 



the expression for the steady-state concentration is universal, i. e. it does not depend on 
solubility laws. 

Let us proceed now to another case: when diffusion fluxes of both dissolved gases are 
significantly non-steady. Here we can use Eq. (I2ip for both i — 1 and i = 2. From Eq. (T27T) 
we have: 

D x -al = D 2 -al (32) 

71 71 

Similarly to Eq. (129|) . we have 

El ( l ~ C l \ 2 = ( ^2,0 ~ 7*2,00 (1 ~ cf) 

D 2 \ c{ ) \ n 1:0 - n 1)(J0 (cf) 

Concentrations and 1 — c\ in the bubble can be found from Eq. f[3"B"j) . when the solubility 
laws [e. g. Eq. (1311 or Eq. (l3"Tj) ] are set. It should be emphasized that the applicability of Eq. 
(1331) is determined only by the condition of validity of self-similar growth of the bubble, i. e. 
the absence of surface tension influence. This condition does not set any limitations related 
with the non-steady character of the diffusion flux of dissolved gas. 

Eq. (I3"3l was not presented in literature before; therefore, the difference between steady- 
state compositions of the bubble for the cases of steady and non-steady diffusion fluxes also 
has not been revealed previously. In Fig. [2] we plot the relative difference between the values of 
steady-state concentration c\ obtained for steady and non-steady diffusion fluxes of dissolved 
Erases as a function of diffusion coefficients ratio D\/D 2 for different supersaturation values. 

In case when both inequalities (1201) and (1211) are violated, equation on the steady-state 
composition can be obtained numerically from Eqs. (|2"T|) . (fill and (13TJ1) or f[3~Tj) . 




3.2 High supersaturations 

If we consider a solution which is strongly supersaturated with both components (such situation 
is typical for homogeneous nucleation of a bubble), the following strong inequality is fulfilled 
for both components 

riifi > n i>tx . (34) 



(35) 



Using inequality fl34j) in Eq. (J29J), for the steady diffusion case we have 

D 2 c{ n lfi ' 

or 



(36) 



i D 2 "i.o ' 
D 1 n 2 ,o 



For the non-steady diffusion case, when Eq. fl33l) is valid, we have 



or, finally, 



(38) 



D 2 "i,o 
Di n 2 ,o 



It has to be emphasized that Eqs. (I3"6"l) and (I3"5|) are universal: they are valid irrespective of the 
solubility law, but only when the solution is strongly supersaturated with both components. 
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4 Conclusions 



While in the case of one-component bubble growth any difference in bubble dynamics for 
different solubility laws vanishes after R exceeds 2<t/IT [14], in the two-component case the 
"memory" of the solubility laws is expressed in the value of the steady-state composition of 
the bubble even when R ^> 2<j/U. 

At the same time, for strong supersaturations the steady-state composition of the bubble 
is universal, i. e. it does not depend on the solubility laws. For strong supersaturations the 
value of the steady-state concentration of each component in the bubble is defined only by 
initial densities of dissolved gases n^ g and diffusion coefficients -D.;. 

Another interesting observation is that the most complex case of bubble growth (when 
supersaturations are extremely high) is the simplest case for multicomponent bubble growth 
conditions. It is the situation when the fluxes are independent from the bubble compositions 
and are determined by the values of components supersaturations. 

Consideration of a two-component gas bubble presented in the current paper evokes the 
following non-trivial question: if the bubble nucleates fluctuationally, its initial (equilibrium) 
composition is defined by thermodynamical parameters [28, 29] and does not depend on dif- 
fusion coefficients Di,D 2 of the components in the solution. But at large bubble radius, as 
was explicitly shown above, the steady-state (but, evidently, non-equilibrium) composition of 
a two-component bubble does depend on these parameters. Therefore, it is important to ob- 
tain the solution to the problem of a physically correct description of a bubble evolution from 
the initial equilibrium composition to the steady-state one. Even for one-component bubbles 
the self-similar approach is not applicable when the bubble radius is small [1, 13, 14] (when 
Laplace forces strongly influence bubble growth). For two-component bubbles the problem is 
even more complex, because the change of composition with time leads to an essential change 
of time dependence of the bubble radius as compared with the self-similar evolution case. 
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0.5 H 




^l,o/^l,oo 

Figure 1: Relative difference between the values of steady-state concentration c\ for Henry's law 
and Sievert's law as a function of n 1)0 /ni jOO for various values of n 2i o/^2,oo : curve I - n 2j0 /n 2iOO = 
2, curve II - n 2>0/ M 2 ,oo = 10, curve III - n 2i o/^2,oo = 100 and curve IV - n 2i0 /n 2>oo = 1000. 
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Figure 2: Relative difference between the values of steady-state concentration c\ obtained for 
steady and non-steady diffusion fluxes of dissolved gases as a function of diffusion coefficients 
ratio Di/D 2 for various values of nifi/ni >ao = n 2>0/ M 2)OO : curve I - ^1,0/^1,00 = ^2,0/^2,00 = 2, 
curve II - ni )0 /n 1)OO = n 2i0/ M 2iOO = 10 and curve III - ni >0 /n 1:OO = n 2 ,o/n 2 ,oc = 100. 
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